A NEW METHOD OF FINDING ALL ROOTS OF SIMPLE 
QUATERNIONIC POLYNOMIALS 



LIANGGUI FENG AND KAIMING ZHAO 



Abstract. In this paper, we provide a new method to find all zeros of 
polynomials with quaternionic coefficients located on only one side of 
the powers of the variable (these polynomials are called simple polyno- 
mials). This method is much more efficient and much simpler than the 
known one in We recover several known results, and deduce several 
interesting consequences concerning solving equations with all real coef- 
ficients or complex coefficients which do not seem to be deduced easily 
from the results in [9]. We also give a necessary and sufficient condition 
for a simple quaternionic polynomials to have finitely many solutions 
(only isolated solutions). 
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1. Introduction 

The quaternion algebra plays an important role in many subjects, such as, 
quaternionic quantum mechanics [T] , and signal precessing [21, IS] . Because of 
the noncommutativity of quaternionic multiplication, solving a quaternionic 
equation of degree n becomes a challenging problem [3-19]. 

Niven in [14:\ [T5] made first steps in generalizing the fundamental theorem 
of algebra onto quaternionic situation which led to the article by Eilenberg 
and Niven [3] where the existence of roots for a quaternion equation of degree 
n was proved using strongly topological methods. After that, Topuridze 
in [20], also with help of topological method, showed that the zero set of 
polynomials with quaternionic coefficients located on only one side of the 
power of the single variable (these polynomials are called simple polynomials) 
consists of a finite number of points and Euclidean spheres of corresponding 
dimension. 

Concerning about the computation of roots of a quaternionic polyno- 
mial, the first numerically working algorithm to find a root was presented 
by Serodio, Pereira, and Vitoria [IS], and further contributions were made 
by Serodio and Siu [TU], Pumpliin and Walcher [T7], De Leo, Ducati, and 
Leonardi [T3], Gentili and Struppa [B], Gentili, Struppa, and Vlacci [^, Gen- 
tili and Stoppato [5j. A large bibliography on quaternions in general was 
given by Gsponer and Hurni in 2006 [Ij. Recently, Janovskd and Opfer 
presented a method in [9] for producing all zeros of a simple quaternionic 
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polynomial, by using the real companion polynomial introduced for the first 
time by Niven [H] , the number one introduced in [16] , and the presentation 
of the powers of a quaternion as a real, linear combination of the quaternion. 
Let us recall some notions first. 

Throughout this paper, let N be the set of positive integers, R the real 
number field, C = M © Ri the complex number field, and H the skew-field of 
real quaternions, that is, any element of H is of the form q = oq + oii + a2i + 
ask = (ao + oii) + (02 + a3i)j, where i, j, k are usual quaternionic imaginary 
units, and ao,ai,a2;fl3 G R; and the R-bilinear product is determined by 

= = k^ = —1, ij = — ji = k, jk = -kj = i, ki = — ik = j. For the 
above quaternion (7, we denote by Reg the real part of g, by \q\ the module 
of q (i.e. \q\ = -^Oq + a| + + Og), and by [q\ the conjugate class of q{i.e. 
[q] = {aqa~^\a G H, o 7^ 0}). We confirm that a quaternionic polynomial 
with the coefficients on the same side of the power of the single variable are 
called a simple quaternionic polynomial or a simple polynomial. 

Definition 1.1. Let zq G M he a zero of a simple polynomial p{z) as 
given on the left-hand side in (2.1) (i.e., p{zq) = 0). If zq is not real and 
has the property that p{z) = for all z £ [zq], then we will say that zq is a 
spherical zero of p{z). If zq is real or is not a spherical zero, it is called an 
isolated zero of p{z). 

In the present paper, we provide a new method for finding all zeros of 
simple polynomials p{z) of arbitrary degree n (Theorems 1 and 4). Our proof 
is based on two well-known techniques: the presentation of a quaternion as 
a 2 X 2 complex matrix, and the Jordan canonical form of a complex matrix. 
We first write p{z) so that its constant term is 1 or 0. Then introduce derived 
polynomials fi{t) and f2{t) oi p{z) which have complex coefficients, where t 
is a real variable, such that p{t) = fi{t) + /2(t)j; and define the discriminant 
polynomial p[t) = fi{t)fi{t) + /2(i)/2(^) oi p{z) which is a polynomial with 
real coefficients, where t is considered as a complex variable, the polynomials 
fi{t) and f2{t) are obtained by only taking the conjugate coefficients of fi{t) 
and f2{t) respectively. Then all zeros ol p{z) can be obtained from complex 
zeros of the discriminant polynomial p{t). More precisely, let G C such 
that p{zq) = 0. If Zq is real then it is an isolated zero of p{z). If zq is not 
real and fi{zQ) = 72(2:0) = /i(-Zo) = f2{zo) = then it is a spherical zero of 
p{z). If Zq is not real and at least one of fi{zQ), f2{zo), fi{zQ), f2{zQ) is not 

zero, let ( ? | be a unit complex solution of the linear system 



Then corresponding to this zq we have an isolated solution for p{t): 



From the above three cases we obtain all zeros of ^(2:). 

The paper is organized as follows. In Section 2, we prove that our meth- 
ods for solving a simple polynomial equation in Theorems 1 and 4 are valid. 
Then we recover several known results, and deduce several very interesting 
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consequences concerning solving equations with all real coefficients or com- 
plex coefficients which do not seem to be deduced easily from the results in 
[9] (see Corollary [2] and Corollary [3]) . We also give a necessary and sufficient 
condition for a simple quaternionic polynomials to have finitely many solu- 
tions (only isolated solutions). In Section 3 we give an algorithm to find all 
zeros of a simple quaternionic equation, based upon our Theorems 1 and 4. 
In Section 4 we give three examples. In particular, we use our method to 
redo Example 3.8 in [9j. In Section 5, we make some necessary numerical 
considerations and compare our algorithm with that in [9]. 

2. Finding all zeros of a simple polynomial 
We consider the simple quaternionic polynomial equation: 

(2.1) g„x" + --- + gix + go = (g„/0), 

where x G HI is the variable, n G N and qi {i = 0, . . . ,n) G H are given. If 
qo / 0, then Eq. (|2.ip can be written as 

qQ^qnX^ H h q^^qix +1 = 0. 

Hence, in order to solve Eq. ()2.ip . it suffices to solve the following equation 

(2.2) p„a;" + ---+pix + (io = 0, 

where pi{i = l,...,n) G IHI,p„ / 0, do = or 1. We simply denote the 
left-hand side of (2.2) as p{x). 

Let o- : M ^ C^^^^ q = zi + ^ { ^— — ] , where zi,Z2 G C. 

\ —Z2 Zl J 

Then a is an M-algebra monomorphism from H to C^^^. Sometimes, this 
monomorphism is also named as the derived mapping of H, and (T{q) is 

a 
a 



denoted by q'^ . Obviously, 



for any a G 



Let Pi 




for i 



1, 



, n, where t^^\t2^ G 



Then 



\ '-2 ''I 

(|2.2p becomes the following matrix equation in the matrix variable Y 



Zl 
-Z2 



Z2 
Zl 



{zi,Z2 G C): 



(2.3) 



t 



(n) An) 



-t 



_ i_ 

(n) ^(n) 
1 



yn + . . . + 



t 



(1) ,(1) 



_ ^ 
(1) ,(1) 



Y + 



do 






do 



0. 



Now we introduce a matrix polynomial P{t) in real variable t (considered 
real number) as follows: 



P(t) (2.4) 



P{t) =pie + ...+p1t + doI, 



where I is the 2x2 identity matrix. 

Write P{t) = Q{t){tl -Y) + Pi{Y), where 

Pi{Y)=p1Y^ + ---+p1Y + doI, 

Q{t) = 
UPiiY) 



,+p1Y)r~' + . . . + (p- +p-y + . . .+p-Y 



n-l\ 



0, then P{t) = Q{t){tl - Y), and hence 

det P{t) = det Q{t) det(t/ -Y) = XY{t) det Q(t), 
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where XY^t) is the characteristic polynomial of Y. 

Set p{t) = detP(t). By Cayley-Hamilton Theorem, we see that p{Y) = 
for every Y satisfying Pi{Y) = 0. 

Notice that, p{t) = det P{t) is a polynomial in real variable t of degree 2n 
with real coefficients, since p{t) = det(p„t" + • • • +pit + doy, and det{pnt^ + 
• • • + Pit + do)'^ ^ for any real value of t. Then 

TU)] (2.5) p^t) = b{t-^if'-^---{t-Csf^-^- 



•(t - - • • • - - 

where 6 is a real number, ^i, . . . are distinct real numbers, r/i,?7i, . . . , 
r]k,f]k are distinct nonreal complex numbers, and ri, r2, • • • , r^, si, • • • , G 
N. It is clear that s + k < n. 

Now suppose y is a solution of Equation (|2.3p . Then, p{Y) = 0. Since y is 

offormf ^ ] where Zl, Z2 G C, then xvit) = {t—zi){t—zi)+Z2Z2 > 

\ ^2 Zl J 

for all real values of t. Consequently, Y has two equal real eigenvalues or two 
conjugate complex eigenvalues. Hence, its Jordan canonical form, Jy, has 

to be of form al or ( 'l ^ ] where a G M and c is a complex with nonzero 



^ c ^ 

imaginary part. Using ()2.5p we know that, for a solution Y of Equation 
(|2.3p . Jy has to be one of the following s + k matrices (up to the order of 
the diagonal entries): 

X3 (2.6) 6/....,«,/,(: ° ).....('? ° 

Next we will prove that each of the above cases can occur. 
Case 1: Jy = ^j/ for i = 1, s. 

In this case Y = Jy- Clearly, y is a solution of Equation (12. 3p iff PnCi" + 
• • • + p^S^i + dol = 0, iff is a common real root of both 



and 



t^a;- + . . . + 41)3; = 0. 



But is a real root of p{t), the equalities 4"'*?? + • • • + 4^'*^* + do = and 
4"'*Cf + • • • + 4^^?« ~ hold naturally. Hence, ^^1 (i = 1, . . . , s) is a solution 



of the equation (|2.3p . 

Case 2: Jy= (^1^ ^ ^ for i = 1, . . . , A;. 

We may assume that (see [22]) 

Y = I ^1 Z2 \ f in \ / Zl Z2 
-Z^Z^jyO Wi j \ -Z2 Zl 



SIMPLE QUATERNIONIC POLYNOMIALS 



eq:2.4 
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for some zi,Z2 G 
and only if 



Pr. 



with l^il^ + |z2| / 0. Then F is a solution of (jO]) if 



Zl Z2 

-z^ zT 











Z2 

-Z^ Zl 
Zl Z2 
Z2 Zl 



-1 



+ 

-1 



I.e., 

(2.7) 



Pn 



Zl Z2 
-Z2 Zl 



m 



+ 



Z2 



+Pl[ - 







- +4 :f 



In other words, for any nonzero solution 
solution for (j2.3p of the following form: 



Wi 

We write ()2.7p into two equations 



Zl Z2 
-Z2 'Zl 









Wi 

Zl Z2 
-'Z2 'zl 

Zl Z2 
-'Z2 Zi 



Zl Z2 
-'Z2 'zl 



0. 



of (|2.7p . we obtain a 



-1 



t^^\^zi + tt'rfli-^) + • • • + ti'^r/izi + t^^^r?,(-z^) + d^zi = 
4''%"^2 + 4" + • • • + 4^%-22 + 4^ W + c^022 = 0, 



I.e., 



eq:2.6 (2.9) -^T 



(4" V + • • • + 4'%)^ + (4"%" + • • • + 4'% + dQ)z2 = 



(4"%" + • • • + 4'% + ^^o)^ - (4"%" + • • • + 4'%)^2 = 0. 



Considering (j2.9p as a linear system in the variables X 
that its determinant of the coefficient matrix is 



Zl 
Z2 



we see 



D 



t\ 'iji"- H r]i 



4"%" + 



+ t'^^^Vi + do 

(trv + ---+4%) 



4"V + 



Since and r/j are roots of 
p{t) = detP(t) : 



t[^)t^ + . . . + 41)^ + do t^^'h^ + - - - ' 



can be given clearly, which will be denoted by Fj. 



t 



+ . . . + 41)^) 4")^" + . . . + + do 

it follows that D = 0, which shows that (|2.9p always has nonzero solutions. 
Since all coefficients of (|2.9p are known, the solution set of (|2.9p in X = 
'zl ^ 

Z2 



lir]i and rji simultaneously satisfy t^"^x"H l-t^^^'x + do 



■■■+t 



(1). 



0, then (2.9) becomes trivial, and Fj = C 



2x1 



Oand 4"^a;" + 
Consequently, 



any element of the form (2.8) is a solution of (2.3). 

Now suppose (2.9) is nontrivial. Then Pj is of dimension 1, and Tj 



6W 



z S C > , in which 



is a fixed nonzero solution of (|2.9p 
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f ormulal 



zerol 



zero2 



with |aW|2 + |&»|2 = 1. Up to now, we have actuahy provided a method to 
find all the roots of (2.1) in H. To summarize our result as a theorem, we 
need to introduce some notions. 



Let Pi = tf^ + G H for i 



1, . . . , n where ^ 



We caU the 



+ t['h + do, f2{t) = tPr + 



following four polynomials the derived polynomials of Equation (2.2) 

m 

flit) - n 

We define the discriminant polynomial of Equation (2.2) asp{t) = /i(i)/i(i)+ 
/2(i)/2(t). We factor it as in (2.5). Remark that these two p{t) are essen- 
tially equal. Introduce sets Ti and T2 as follows 

Ti = {77 e {vi, ...,%} I fiiv) = f2{ri) = Mv) = Mv) = 0}, 

T2 = {m,---,Vk}\Ti. 

Now we can state our main result 
Theorem 1. With the above notations, the solution set of (2.2) overM. is 



(2.10) 

where oji takes 
(2.11) 



{^1 , . . . , 6}U,,,eT2 {wi} U^^eTi h 



1 



:{\f2ir]i)\% + \fl{r]i)\'rn-2hiv^)fliv^)il^m)k} 



\fliV^W + \f2{r^^)\' 

as its value |/i(??j)P + l/2(^i)P / 0, otherwise takes 

(2.12) ^— r^^^— r^{|/i(r^)|N. + |/2(r^)P???+2/2(W)^M(Imr/0k} 

l/i(^i)l + \j2{vi)r 

as its value, in which Imrji means the imaginary part (real number) of r]i. 
Moreover, the union of the first two parts in the above set is the set of isolated 
solutions and the third part of the above set is the set of spherical solutions. 

Proof. To finish the proof we need to continue the argument on the case 

6» ) 
X = with 



when (2.9) is nontrivial, i.e., rji £ T2. For any rji £ T2, let 
a nonzero complex solution of the system 



f2{vi) him. 
flini) -f2iv 

Since the set of nonzero solutions of (j2.9p in unknown X 



Zl 
Z2 



is <^ / 



/ 7^ 0, Z e C > , the solutions of (|2.3|) corresponding to 



the Jordan canonical form 



Vi 








are 




n(0 
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which is actually one value, and in H which can be written as 

which can not be a real number. Therefore, the solution set of Pnx"^ + 
pix + do = in the skew- field H is 
(2.13) _ 

{6, . . • ,6}U^,eT2 {W^'^W + |6«|2r^- 2a«6»(Imr?,)k} U,,^eTjr?i]. 

Note that |/i(r/i)|2 + |/2(7?i)|2 and + |/2(t^)P can not be simulta 

neously for rji G T2, then we can take 



+ 



6« 



Vl/i(^)P+l/2{^?r)P 

V Vl/l(W)P + l/2{W)P 



for the case + |/2(??j)l^ 7^ 0) otherwise we take 



hiM) 



^\h{m)?+\f2{mW 
__hm__ 



After manipulations, the set in (|2.13p becomes 
where a;,- takes 



1 



-{|/2(^?.)r^?. + |/i(r?i)|"fy-- 2/2(r/i)/i(r/i)(Im??.)k} 



l/i(^.)P + l/2(r/.)|2 
as its value if + l/2('?j)P 7^ 0; otherwise it takes 



1 



7=7^{|/i(W)r^^ + l/2(r?l)rr?-+ 2/2(r^)/i(r^)(Im,?,)k} 



|/i(W)P + l/2(r^)P 
as its value. 

Finally, it is clear that [r/j] contains no real numbers for r/,, € Ti. This 
completes the proof. □ 

Note that there is no repetition in the solution set given in (j2.10p and one 
can use (121111 or (l2T2]l for uji if 

mm)? + \f2{m)?){\h{m)\'' + \f2{m)?) 

Theorem [1] shows, once we get a decomposition with the form ()2.5p of the 
discriminant polynomial then we can produce all roots of the quater- 
nionic equation (|2.2p by formula ()2.10p . 

From Theorem 1 and the arguments before Theorem 1, we can easily see 
the following well-known results. 



corollary2 Corollary 2. (a). Any simple quaternionic equation QnX^ -|- • • • -|- qix + 
go = 0(g„ 7^ 0) has a root in H. 
(h). The simple quaternionic equation qnX^ -!-•••+ qix + qo = 0(g„ 7^ 0) 
has a finite number of roots in EI iff it has at most n distinct roots 
in m. 
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(c). The roots of QnX^ + • • • + qix + qo = 0(g„ 7^ 0) are distributed in at 
most n conjugate classes, and there are at most n real roots among 
them. 

Proof, (a) is obvious. 

(b) From Theorem 1 we need to only show that [r]i\ is an infinite set if 
r]i G Ti. It is weU-known from [21] that mi,U2 £ H are conjugate (i.e., 
there exists nonzero g G HI such that ui = qu2q~^) iff Reui = Re'U2 and 
\ui\ = 1^2 1 . Since rji £ Ti has a nonzero imaginary part, [rji] is an infinite 
set. Thus qnx"^ + • • • + qix + qo = has a finite number of roots in EI iff 
Ti = 0, iff PnX^ + • • • +pix + do = has at most n roots in EI since s + k < n 
(See ()2.5p for the notations). 

(c) fohows from s + k <n. □ 

Now we give a quick method to solve simple quaternionic polynomials 
with all real coefficients or with all complex coefficients. These results do 
not seem to be deduced easily from the results in [9j. 



corollaryS Corollary 3. (a). If all qi in qnx"^ + • • • + qix + qo = 0{qn ^ 0) are real 
numbers and the solution set of this equation in C is {^1, . . . , Cs; Ci) Ci) 
. . . ,Q,Ct}, where ^1, . . . , o'^e distinct real numbers, Ci) • • • 1 Ct 0,''"^ 
distinct nonreal complex numbers, then the solution set of this equa- 
tion in M is 

{ei,...,es,}u[Ci]u---u[Ct]. 

(b) . More generally, if all qi in qnX^ + ••• + qix + qo = 0{qn 7^ 0) 

are complex numbers, and the solution set of this equation in C is 
{Ci,---,^s,Vii---,mXi,Ci,---,Ct,Ct}, where ^i,...,^s are distinct 
real numbers, rji, . . . ,rjk,Ci, ■ ■ ■ ,Ct are distinct nonreal complex num- 
bers (eachrji is no longer the root of this equation), then the solution 
set of this equation in EI is 

{6,...,6,m,---,%,}u[Ci]u---u[Ci]. 

(c) . Let fi{t), f2{t), flit) and f2{t) be the derived polynomials for 112. gj) . 

Then \2.2) has finitely many solutions iff the complex polynomial 
gcd(/i(t), /2(t), /i(t), /2(t)) has no nonreal complex root iff the com- 
plex polynomial gcd(/i(t), /2(t)) has no nonreal conjugate complex 
roots. 

Proof, (a) This is a special case of Part (b). 

(b) When the equation considered has only complex coefficients, two of 
the derived polynomials are fi{t) = qnt"^ + • • • + qit + qo (up to a com- 
plex scalar), and f2{t) = 0. So, the roots of the discriminant polynomial 
are ^1, . . . , m> • • • , ^fc,^, Ci, Ci, • • • , Ct, G- For each r]i{i = l,...,k), 

( a(*) \ / 1 \ 
since f \(r\i) / 0, is in T2. We can take ( ^j^^-j j = ( ^ j as the 

unit complex solution of the system { {^J^ W^ V = 0. Then, 

|aW|2?7j + |6W|27^-26WaW(Imr/i)k = It is easy to see that e T^. This 
completes the proof. 
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(c) Suppose c is a nonreal complex root of gcd(/i(t), /2(t), /2(t)). 
Then both c and c are roots of the discriminant polynomial, /i(c) = /2(c) = 
flic) = /2(c) = and /i(c) = /2(c) = /i(c) = /2(c) = 0, which imphes_that 
at least one of c, c is in Ti. Thus, Ti = iff gcd(/i(f), /2(t), /2(t)) 
has no nonreal complex root. From Theorem 1 we see that ()2.2p has finitely 
many solutions iff Ti = 0. The conclusions in the corollary follow easily. □ 

Now we can give a simplified version of Theorem 1. 

Theorem 4. Consider the simple quaternionic equation p{x) := PnX^ + 
• • • + pix + do = 0, where € IHI with Pn ^ and do = or 1. We 
write pit) = git)igiit) + 52(i)j) where t is considered as a real variable, 
9i9ii92 £ with gcd(5i,g2) = 1- Let the complex solution sets for git) 
and git) = + g2(i)52(t) are 

(2.14) {6,---,6;Ai,Ai,A2,A2,--- ,At,At;r/i,...,?7fcJ, 

(2.15) {r/fc^+i,r?fc^+i . . . ,r?fc,r?fc} 

respectively. We may assume that ^1 , . . . , are distinct real numbers; rji, . . . , 
%i,%i+i,%i+i---,%,%,Ai,Ai,A2,A2,--- ,At,At are d istinct nonreal com- 
plex numbers after deleting those rji 's and fii 's in 12.15\) if one of them ap- 
pears in the set ^2.14^ . Then the quaternionic solution set of pix) is 

(2.16) 

where 



bi('7i)P + IS'2(??i)P 



{\92iVi)\ Vi + \9iiVi)\ r]i + '^92ir]i)giiVi)i^T0tir]i)k}. 



Proof. We continue to use the notations in Theorem 1. We know that /i(t) = 
9it)9iit), f2it) = 9it)92it), Pit) = git)git)git), and git) has no real root. 

From Theorem 1 we see that {.^1, . . . , ^s}UjLi[Ai] are zeros of p(x). All 
other zeros come from {r/i, . . . , 7/^^ , r^^i+i, %i+i • • ■ ,??fc, ??fc}- For each iji, we 
see that |/i(??j)P + |/2(f?i)P 7^ 0. Then using (|2.12p in Theorem 1 and 
simplifying we obtain Wj. This completes the proof. □ 

Remark that the above theorem simplifies the computation for finding all 
zeros of p(x), and the following known result (see [I6J) can follow easily from 
the above theorem. 

Corollary 5. The spherical zeros of simple quaternionic polynomial pix) := 
PnX^'' + • • • + pix + poipi G M,pn ^ 0) are distributed in at most INT(^) 
conjugate classes, where INT(^) means the integral function value at ^. 

3. Algorithm 

Based on our Theorem 1, we now can give an algorithm to solve the 
quaternionic equation qnX^ + ■ ■ ■ + qix + qo = iqn / 0), as follows. 



Algorithm 1 (for solving the simple quaternionic polynomial equation 



H h gix + 50 = 0) 

Step 1. Write the equation as Pnx'^ + • • • + pix + do 
or 1 (in fact, if qo ^ 0, simply multiply the equation by 



with do = 
on the left). 
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Write Pi = t^l^ + 4*^j for i = 1, 
polynomials and discriminant polynomial of Pnx'^ + • • • + pix + do = 0: 



, n with t^\t2^ G C. Find the derived 



= + . . . + t'^'t + do, f2{t) = t^i" + • • • + ^2''*; 



f,{t) = t^^'h^ + ■■■ + t\'>t + do, f2{t) = t^^'f" + ■■■ + A^'t, 

P{t) = fl{t)m+f2{t)f2{t). 

Make sure the coefficients of p{t) are real. 

Step 2. Compute all distinct zeros (real or complex) of the discriminant 
polynomial p{t) (in MATLAB, use the command roots). Denote these ze- 
ros by ^1, . . . , ^s, rji, . . ., rjk, rji, . . ., rjf. such that ^1, . . . , are distinct real 
numbers and rji, . . ., r]k are distinct nonreal complex numbers. Then classify 
rji, . . . ,rik into two sets Ti and T2 : 

Ti = {r/ G {7?i, ...,%} I /i(r?) = ^(r/) = hiv) = Uv) = 0}, 

T2 = {m, - ■ ■ ,Vk}\Ti. 

When the simple quaternionic polynomial considered is a polynomial with 
real coefficients, then directly set Ti = {r]i, . . . , r/^} and T2 = 0. 

Step 3. For each r]i £ T2, compute fi{rji), /2(??i)- Then by Formula 
(|2.10p . output all roots of g'nx" + • • • + qix + go = 0. 

If we use our Theorem 4, then we get a better version of Algorithm 1. 



Algorithm 1 ' (for solving the simple quaternionic polynomial equation 
Qnx" H h gix + 50 = 0) 

Step 1. Write the equation as p{x) := PnX^ + • • • + pix + do = with 
do = or 1 (in fact, if qo 7^ 0, simply multiply the equation by q^^ on 
the left). Then write p{t) = g{t){gi{t) + g2{t)i) where t is considered as 
a real variable, g,gi,g2 £ C[i] with gcd{gi,g2) = 1. Now we compute 
m = 9i{t)gi{t) + g2{t)g2{t). 

Step 2. Compute all distinct zeros (real or complex) for g(t) and g{t) 
respectively (in MATLAB, use the command roots): 



first2 (3.1) 


{6) ■ ■ ■ ,^s; Ai, Ai, A2, A2, • • 


■ ,At,At;r/i,. 




sec2 (3.2) 


{Vki+l, Vki + 1 ■ ■ 







where ^i,...,^^ are distinct real numbers; {Ai, Ai, A2, A2, • • • ,Xt,\t} and 
{r]i, . . . ,r]ki, r]ki+i,f]ki+i • • • , i]k,Vk} are two sets of distinct nonreal complex 
numbers. Delete those r/j's and fji's in (|3.2|) if one of them appears in the 
set (IXTTl . 



f ormulaomega 



f ormula21 



Step 3. For each r]i, we compute 
(3.3) 



1 



:{\92{Vi)\ Vi + \9i{Vi)\ Vi + '^92{Vi)9iiVi)0-T^Vi)^- 



\91iViW + \92iViW 
Then the quaternionic solution set of p{x) is 
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Prior to our method, Janovska and Opfer gave an algorithm in [5] for 
solving the same quaternionis equation QnX"^ + • • • + qix + qo = {qn 7^ 0). 
In order to introduce their algorithm precisely, let us first recall the known 
concept, companion polynomial. For 



compciaion 



Pix) = ^qjx^,qj e M,j = O, . . . ,n, go,9n 7^ 0, 

j=0 

following Niven |14j , or more recently Janovskd and Opfer [9j , its companion 

polynomial is defined by 

(3.5) 

n 2n min(k,n) 

q2n{x) = ^ qjqkX^^^ = y^fcfca:^, where hk = ^ ^jlk-j G 
j,k=0 k=0 j=max{0,k—n) 

We remark that q2n is equal to the discriminant p{t) in our case. From 
Pogorui and Shapiro|16j, we know that all powers x^ , j G Z of a quaternion x 
have the form = ax+P with real a, (3. In particular, x'^ = 2(Rex) x — 
In order to determine the numbers a, /3, Janovskd and Opfer in [9] set up 
the following iteration: 



eq: iterationl (3.6) 



x^ = ajX + f3j , , /3j E M, j = 0, 1, . . . , 
ao = 0,/3o = 1, 
Oj+i = 2Rex Oj + f3j, 
/3j+i = -kpoj, j = 0,1, ... . 



Now by means of the first line of Iteration ()3.6p . the polynomial p{x) can be 
rewritten as 



formulas 



P{^) = "^QjioijX + /3j) = "^qjaj \ x + I "^Qjf^j ) = A{x)x + B{x), 
3=0 \i=o / \j=o 

where 

n n 

(3.7) A{x) = ^qjaj, B{x) = ^nj^j. 

j=0 j=0 

With these in hand, we can state the algorithm given in [5], as follows. 

Algorithm 2 (for solving the simple quaternionic equation qnx'^ + • • • + 
qix + qo = 0) 

Step 1. Write + - • ■+qix + qo = asp(x) = anX^ + - ■ ■+aix + aQ = 
with a„ = 1. For this p{x), compute the real coefficients bo, bi, . . ., b2n of 
the companion polynomial q2n{x) by formula (|3.5p . 

Step 2. Compute all 2n (real and complex) zeros of q2n{x), denote these 
zeros by zi, Z2, • • •, Z2n and order them (if necessary) such that -Z2j-i = ztJ, 
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j = 1,2, .. ., n. 

Step 3. Define an integer vector ind (like indicator) of length n, and 
set all components to zero. Define a quaternionic vector Z of length n, and 
set all components to zero. 

For j := 1 : n do 

(a) Put z := Z2j-i. 

(b) if z is real, Z{j) := z; go to the next step; end if 

(c) Compute v := A(z)B{z) by formula ()3.7p . with the help of Iteration 
HIS]). 

(d) if w = 0, put ind(j) := 1; Z{j) := z; go to the next step; end if 

(e) if f 7^ 0, let {vi,V2,V3,V4) := v. Compute \w\ := -^^2+^3+^4, and 
put 



f ormula4 



(3.8) 

end if 
end for 



ZU) := Rez 



llmzl 



llmzl 



llmzl 



\w\ 



-V2,- 



W 



-V3, 



W 



In this algorithm, corresponding to a real z the expression Z{j) produces a 
real isolated zero z, corresponding to "f = 0" the expression Z(j) produces 
a spherical zero [z], and corresponding to "f 7^ 0" the expression Z{j) 
produces an isolated zero Rez — ■'^^^'21 — -^^^p^^sj — "^^^ output 

results of Z{j) produce all zeros of polynomial qnx'^+- ■ ■+qix+qQ. Algorithm 
2's original edition is Section 7 of [9j, where \w\ := \/v\+ v'^ + v\ is false, 
which is a misprint. 

4. Examples 

Example 1 In H, solve the equation p{x) := ix^ + jx^ + kx + 1 = 0. 

Solution 1. Use our Algorithm 1 to do this. Write the coefficients i,j,k 
into the form: i = i+Oj, j = 0+lj, k = 0+ij. Then the derived polynomials 
of this equation are 

/l(t)=it3 + l, /2(t)=t2 + it, 
/l(t) = -it3 + l, Ut)=t^-\t, 

the discriminant polynomial is 

pit) = hh + f2h = {t-i){t + i)(t - e'5)(t - e'x)(i - e'x)(t - e'f), 

which has no real root. It is easy to see that 

Ti = and T2 = {rji = i, 7^2 = e * , 7^3 = e * }. 

For r]i = i, then /i(i) = 2 / and /2(i) = —2, and we get an isolated 
zero by Formula (|2:T0D : | • (4i - 4i - 2 • (-2) • 2 • 1 • k) = k; 

For r?2 = e^f = ^ + ^i, then h{r]2) = ^ - ^i and ^(r/a) = + 

"^"^^ i, and we get another isolated zero by Formula (j2.1U|) : + |i + |k; 

Similarly for 773 = we get the isolated zero: + ^i + ^k. 
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Thus, the solution set is |k, ^ + + ^k, + + ^k|. 

Solution 2. Use our Algorithm 2 to do this. Write = (it^ + 1) + if' + 
:= /i + /2j, then 

9 = gcd(/i, f2) = t + i, c/i = + t _ i, ^2 = 5 = 51^1 + 5252 = t^ + 1, 

Compute the zeros of g: r/i = — i. 

Compute the zeros of ^: ?/2 = e'*, r/2, % = e'~, 7^3. 

Now for ?7j(i = 1,2,3), compute |gi(r/l)P, 5'2(??l), and |52(??i)P, by 

formula ()3.3p we get the solution set of ix'^ + j^;^ + kx + 1 = 0: 

f, ^/2 1. 1 ^2 1. 1, 1 
< k, — + -i + -k, -— + -i + -k ^ . 
1'2 2 2'2 2 2j 

Example 2 Solve the equation x'^ + + x + 1 = in H. 

Solution 1. We will use our method to do this first. Since x^+x^+x+1 = 
(x + l)(x + i)(x — i), from Corollary [3ja) , we directly know the solution set 
is {-l}U[i]. 

Solution 2. Now we use the method in [9] to do this. 
First Step. By formula (|3.5p . we get 

qQ{x) = (x^ + x^ + X + Vf = + 2x^ + 3x'' + 4x^ + 3x2 + 2x + 1. 
Second Step. Compute all zeros of (2;): — 1, — 1, i, i, — i, — i. 
Third Step. For —1 we get a real isolated zero —1. For i we need give the 
following expansion by Iteration ()3.6p : 

i° = Oi + 1, = 11 + 0, = Oi + (-1), = (-l)i + 0. 

Then by formula ()3.7p get 

^(i) = + 1 + + (-1) = 0, B{i) = 1 + + (-1) + = 0, 

and V := A{\)B{i) = 0. It leads to a spherical zero [i]. Finally, for — i we 
have to repeat the same process as that made for i: 

(-i)0 = 0(-i) + l,(-i)i = l(-i) + 0, 

(-i)2 = o(-i) + (-l),(-i)3 = -l(-i)+0, 

^(-i) = + 1 + + (-1) = 0, B{-\) = 1 + + (-1) + = 0, 

and V := A[—\)B{—\) = 0. It also produces a spherical zero [— i]. Note that 
[— i] = [i], so the solution set is {— 1} U [i]. 

At last let us solve the same polynomial in Example 3.8 of 

Example 3. Find all zeros oi p{z) = + j-z^ + iz^ — — jz — i. 

Solution 1. We first use the method in this paper to do this. 
We have ip{z) = iz^ + ijz^ — — iz^ — ijz + 1. Then 

ipit) = (it^ -t^- it^ + 1) + (it^ - it)3, 

9 = - 1, 91 = - 1, 92 = it, 
9 = 9191 + 9292 = {it^ - l)(-it^ - 1) + it • (-it) = t^ + t^ + l. 



14 LIANGGUI FENG AND KAIMING ZHAO 

The all distinct zeros (real or complex) for g and g are respectively {1,-1, 
and {e~'t,e't,e~'T",e'"f }. Now we take 



rji = e-'3 = 1/2 - 
rj2 = e~'^ = -1/2 - 

giim) = ie'^ - 1 = -1 - ^^ - li, |gi(f/i)|2 = 2 + ^/3; 

52(fyi) = ie'f = -^ + ii, |52(j^i)|2 = 1; 

Sife) = ie'^ - 1 = -1 + ^ - ii, Igiim)? = 2-^/3; 

52(r?2) = ie'^ = - ii, b2(??i)P = l. 
So by formula p.3p . we get 



W2 



: 1 1 If 

2 2 2"' 2 



' K(4 + ^i) + (2-^/3)(-l-^i) 



3-V3 2 2 ' ' 2 2 

+ 2(-^-ii)(-l + ^ + ii)(-^)k} 
^2 2 2 2 2 ^ 

1 1. 1. 1, 

= h -1 J k. 

2 2 2"' 2 

Hence, the set of roots of p{z) is 

' ^ ^2 2^ 2"^ 2 ' 2 2^ 2"^ 2 ^'^^ 



Solution 2. We use the method in j9] to redo this. 

Step 1 Compute the companion polynomial 512(2;) by formula (|3.5 

qi2{x) =x^^ +ix^^ + - - jx^ - ix*^ 

- jx" + x^° - jix'^ + jx"^ - x*^ + jix^ 

- ix^° - kx^ + x*^ + ix*" + kx^ - x"' 

- X — JX — IX + X + JX + IX 

+ ix^ — x^ + jix^ — jx^ + x^ — jix 
+ ix^ + kx^ — x^ — ix^ — kx + 1 



=x^2 + + x^° + - x^ + - 2x^ + - x'^ + + x^ + 1 
=x^2 + x^° - x^ - 2x^ - x'' + x^ + 1. 
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Step 2 Compute the zeros of qi2{x). 



15 



(li2{x) =x^{x'^ + — 1) — a 

=x8(x^ + X^-1)- x2(x^ + x2 - 1) - X^ + 1 

={x^ + x^ - l)(xS - x^) - (x^ - 1) 
=(x^-l)(x2 + l)(x6-l). 
The 12 zeros of q\2 are 



11 1 1 • • • . 1 1 1 

1, 1, —1, —1, 1, 1, —1, —1, — 1, 1, 1, 

'' ' ''' ' '222 222 



Step 3 For 1 and —1, we get isolated zeros: 1,-1. 



1 
2 



^/3. 



For X 



i, then 

[0 = Oi + 1 
1 = li + 

[2 = 0i + (-l) 

" = (-l)i + 
[4 = Oi + 1 
.5 = li + 
^6-0i + (-l) 

A{{) = 1-0+j- 1 + i- 
+0-(-l) + (-l)-0 
+(-j)-l + (-i).0 = 



)0 = o(-i) + 1 
)i = l(-i) + 

)2 = 0(-i) + (-l) 
)3 = (-l)(-i)+0 

Y = o(-i) + 1 

)5 = l(-i) + 

)6 = 0(-i) + (-1) 







A(-i) = 1-0+j- 1 + i 
+0-(-l) + (-l)-0 
+(-j) - 1 + (-i) -0 = 







V = A{i)B{i) = v = A{-i)B{-\)=0. 

They produce the same spherical zero [i] since [— i] = [i]. 

For X = 5 + ^i, 5 - ^i, we have 



(^+ 


2 






) + l 


/I v^, 

V2 2 




= 0(i-fi) + l 


a+ 


V3 
2 






) + o 


^1 ^3 
V2 2 ' 




= l(i-f i) + 


ih+ 


V3 
2 




= l(i + ^i 


) + (-!) 


(1 V3, 
\2 2 ' 


? 


= i) + (-l) 




V3 
2 




= 0(i + ^i 


) + (-!) 


(1 V3 
V2 2 


f 


= 0(i-f i) + (-l) 




V3 
2 




= (-l)(| + 


^i) + 


(1 V3 
*.2 2 ' 


r 


= (-l)(i-^i) + 


(1+ 


%/3 
2 




= (-!)(! + 


fi) + l 


(I \/3, 
V2 2 ' 




= (-l)(i-^i)+l 


(1+ 


2 


i)6 


= 0(i + ^i 


) + l 


(1 Vs. 

^2 2 ' 




= 0(i-^i) + l 



^(a;) = 1-0+j- -1 + i- -1 
+0 - - 1 - 1 + -j - 1 + -i - 

= -l-i-2j 

B{x) = 1- 1+j - 1 + i-O 
-l--l + -j-0 + -i-l + 0--l 
+0--l = 2-i+j 



A{x) = 1-0+j- - 
+0 - + • • • 

= -l-i-2j 

5(x) = l-l+j-l 
+i - + - - - 
=2-i+j 



l + i--l 



V = A{x)B{x) = -3 + 3i + 3j + 3k v = A{x)B{x) = -3 + 3i + 3j + 3k 
\w\ = 3\/3 \w\ = 3\/3. 
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They produce the same isolated zero: | — ■ 3i — ■ 3j — ■ 3k, which 
is 0.5 - 0.51 - 0.5j - 0.5k. 



Finally for 


X — 2 


+ 


2 ^' 


1 

'2 


2 ^' 


we 


still have to finish the followinj 


work: 
















V 2 ' 2 ' 


i)"=0(- 


1 

" 2 ' 




+ 1 




(- 


_1_ ^^i)0 = 0(-i-^i) + l 


(-^ + # 


i)l = l(- 


1 
2 




+ 




i' 


-i-^i)i = l(-i-^i)+0 




= -l(- 






-(- 


1) ( 


1 

2 




V 2 ' 2 ' 


i)3=0(- 


1 
2 




+ 1 




i- 


-i-fi)3 = 0(-i-fi) + l 




ir = i(- 


1 

" 2 ' 


+ ^i) 


+ 




i- 




(-1 + ^1)5 = 




1 

2 " 


+ #i) 


+ (- 


-1) (- 


1 

2 






i)'^=0(- 


1 
2 


+ #i) 


+ 1 




(- 


^^i)6^0(-i-^i) + l 



= l-0+j--l + i- 1 + 0-0 A{x) ==l-0+j--l + i-l 
-1 • -1 + -j • 1 + -i- +0-0 + --- 

= 1 + i - 2 j = 1 + i - 2 j 

B{x) = 1 • 1 + j • -1 + i- + 0- 1 B{x) = 1 • 1 + j • -1 + i - 
-l--l + -j-0 + -i-l +0-1 + --- 

=2-i-j =2-i-j 



V 



A{x)B{x) = 3 - 31 + 3j + 3k v = A{x)B{x) = 3 - 31 + 3j + 3k 



\w\=3VS |w|=3%/3. 
By formula (|3.8p . They produce the same isolated zero: —0.5 + 0.51 — 0.5j 
0.5k. The set of roots of p{z) is 



{1, -1, 0.5 - 0.51 - 0.5j - 0.5k, -0.5 + 0.51 - 0.5j - 0.5k} U [1]. 



In Example 3.8 of [9j, the authors said that the 12 zeros of qi2 were 1 
(twice), -1 (twice), ±1 (twice each), 0.5(=blibi). We point out that 0.5(±libl) 
should be corrected as 0.5(ibl it \/31). 



5. Numerical consideration and algorithm comparison 

The polynomial in Example 1 has the property that its discriminant poly- 
nomialp can be factored by hands as (t— l)(t+l)(t— e'* )(t— e'~)(t— e'~)(t— 
e'~). In general, one can not always expect to get the zeros of the discrim- 
inant polynomial by hands, one has to rely on machine computations. In 
Algorithm 1, if we compute the zeros of the discriminant polynomial, p{x), 
of Example 3 in Section 4 by MATLAB 7.12.0(R2011a), we find the 12 zeros 
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are as following: 



(5.1) 



1 
1 


1 nnnnnnnnnnnnnni 

— l.UUUUUUUUUUUUUUi 


_i_n 


nnnnnnnnonfiRf^/i o^' 

UUUUUUUUZUDD04Zi 


o 

I 


—l.UUUUUUUUUUUUUUi 


— U 


UUUUUUUUzUDD542i 


3 


-U.5UUUUUUUUUUUUUU 


+u 


866U254U378444Uii 


4 


-U.5UUUUUUUUUUUUUU 


-u 


866U254U378444Ui 


5 


U.99999999U1U23U4 


+u 


UUUUUUUUUUUUUUUi 


6 


1.UUUUUUUU9897694 


-u 


UUUUUUUUUUUUUUUi 


7 


U.5UUUUUUUUUUUUUU 


+u 


866U254U3784439i 


8 


U.5UUUUUUUUUUUUUU 


-u 


866U254U3784439i 


9 


U.UUUUUUUUUU16U75 


+1 


UUUUUUUU8531U5H 


10 


U.UUUUUUUUUU16U75 


-1 


UUUUUUUU8531U5H 


11 


-U.UUUUUUUUUU16U74 


+u 


999999991468949i 


12 


-U.UUUUUUUUUU16U74 


-u 


999999991468949i 



which are nearly same as Table 1 in [9]. This is simply because our p{t) is 
exactly the companion polynomial qi2{t)- Hence, the similar measures to 
[9] can be made to obtain machine precision for the zeros with multiplicity 
2(e.g., by Newton's method). It is interesting to note that the discriminant 
polynomial given in this paper is always equal to the companion polynomial 
after considering the variable as a real variable. 

In Algorithm 1, if we have found by MATLAB the zeros of discriminant 
polynomial, we usually by comparing |/i(f/i)P + |/2(^i)P with |/i(r/i)p + 
|/2(^)P decide which one u)i should take between (|2.1ip and (|2.12p . If 
|/i(^i)P + |/2(??i)P is greater, then LOi takes ()2.1ip . otherwise takes (|2.12p . 
However, our Algorithm 1' does not need to make such a decision. 

In Step 2 of Algorithm 1, we need decide whether a zero ^ is real. In our 
experience (the same as [9J), a test of the form |Im^| < 10^^ is appropriate. 
There is another delicate decision to make in our Algorithm 1. That is, in 
Step 2 one has to decide whether a zero r/ satisfies /i(r/) = /2(^) = fiiv) = 
f2{v) — 0- For this, a test for /(?]) = can be carried out in the form 
|/(ry)| < lU~"^''(This was also used to test f = U in [9j). Our Algorithm 1' 
avoids making this delicate decision. 

Both Algorithm 1 and Algorithm 2 need to find all zeros of the compan- 
ion polynomial(note that the discriminant polynomial is always equal to the 
companion polynomial). Unlike Algorithm 2, our algorithms no longer need 
the use of any iterations. Using our algorithms, one can easily produce all 
zeros of the simple quaternionic polynomial from the zeros of the companion 
polynomial. Our algorithm has less workload, as shown in Section 4. This 
has a great advantage for a simple quaternionic polynomial with high de- 
gree. In fact, when one use Iteration (|3.6p to compute A(z) and B(z) for a 
nonreal complex root z with \z\ ^ 1, if the degree of the simple quaternionic 
polynomial considered is high, the workload is huge, even interferes one's 
deciding whether v is zero. 

In [Page 252, 9] D. Janovskd and G. Opfer said "We made some hundred 
tests with polynomials pn of degree n ^ 50 with random integer coefficients 
in the range [—5, 5] and with real coefficients in the range [U, 1]. In all cases 
we found only (nonreal) isolated zeros z. The test cases showed |pn(^)| ~ 
10~^^. Real zeros and spherical zeros did not show up. If n is too large, 
say n ~ 100, then usually it is not any more possible to find all zeros of the 
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companion polynomial by standard means (say roots in MATLAB) because 
the coefficients of the companion polynomial will be too large." But this can 
be easily avoid by using Algorithm 1' in many cases. For example, to solve 
a real coefficient polynomial of degree of 50, we can use Corollary 3 to do 
this easily. 

Example 4. Find all zeros of p{z) = z^^^^ — 2 in H. 
Solution. We know that the complex solution set of p{t) is 

r I 1000/7^1 , ,499 r looo/?: hjil _ 

{zb V2} U^Li im = V2e5oo , T]^}. 

Applying Corollary 3, we obtain the solution set for p{z): 

{±l}ufii["°^eM]. 

According to the above comments D. Janovskd and G. Opfer made on 
their method, it is impossible to find all solutions of this example since the 
degree of p{z) is much larger than 100 and all roots except for it ^° are 
spherical. 

6. Conclusion 

The method introduced in Section 2 can also be used to find all zeros of 
simple polynomials with quaternionic coefficients located on only right side 
of the power of the variable, since x"g„ + • • • + xqi + go = is equivalent to 
q^af^ + ■ ■ ■ + qix + qo = 0. In retrospect, Eilenberg and Niven in [3] for the 
first time proved that quaternionic qnX^ + • • • + qix + qo = has always a 
root, but the method they used is topological. For this result, here we have 
provided a new alternate proof, without using any topological tools. Based 
on the derived polynomial and discriminant polynomial introduced in this 
paper, we have given the zeros in an explicit form for a simple quaternionic 
polynomial. Comparing with the methods provided in ^ to seek roots of 
qnX^ + • • • + qix + qo = 0, our method is different, and has great advantages. 
Using our method, we have recovered several known results, and deduced 
several very interesting consequences concerning seeking the zeros of a simple 
quaternionic polynomial which do not seem to be deduced easily from the 
results in P]. 

We would like to conclude our paper by making the following remark. We 
hope the method in the proof of our main theorem can be useful to seeking 
a method to find all solutions of a quaternionic polynomial equation with 
coefficients on both sides of the powers of the variable. 
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